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Halfvortices in flat nanomagnets 

Gia-Wei Chernp David ClarkeLlJ Hyun Youk||| and Oleg Tchernyshyov 
Department of Physics and Astronomy, The Johns Hopkins University, Baltimore, Maryland, 21218, USA 

We discuss a new type of topological defect in XY systems where the 0(2) symmetry is broken in 
the presence of a boundary. Of particular interest is the appearance of such defects in nanomagnets 
with a planar geometry. They are manifested as kinks of magnetization along the edge and can be 
viewed as halfvortices with winding numbers ±1/2. We argue that halfvortices play a role equally 
important to that of ordinary vortices in the statics and dynamics of flat nanomagnets. Domain 
walls found in experiments and numerical simulations are composite objects containing two or more 
of these elementary defects. We also discuss a closely related system: the two-dimensional smectic 
liquid crystal films with planar boundary condition. 



I. INTRODUCTION 

It is well known that topological defects play an im- 
portant role in catalyzing the transitions between differ- 
ent ordered states for systems with spontaneously broken 
symmetries [l|,[4|. For instance, in nanorings made of soft 
ferromagnetic material, the switching process usually in- 
volves creation, propagation, and annihilation of domain 
walls with complex internal structure 0, H]- We have 
pointed out in a series of papers p-0] that domain walls 
in nanomagnets of planar geometry are composed of two 
or more elementary defects including ordinary vortices in 
the bulk and fractional vortices confined to the edge. The 
simplest domain wall in a magnetic strip consists of two 
edge defects with opposite winding numbers n — ±1/2. 

In a nanomagnet with the geometry of a disk, the 
strong shape anisotropy due to dipolar interaction forces 
the magnetization vector M to lie in the disk plane, effec- 
tively making the magnet a 2D XY system. At the edge 
of the film, dipolar interaction further aligns the spins 
to either of the two tangential directions of the edge: 
rii = M/|M| = ±T. The reduction of ground-state sym- 
metry from 0(2) to a discrete Z2 allows for a new type of 
topological defect confined to the edge. These edge de- 
fects are manifested as kinks in magnetization m along 
the boundary. Kinks are topological defects connecting 
different types of degenerate ground states in ID sys- 
tems with discrete symmetries such Ising ferromagnet; 
their topological properties usually are rather simple [l|. 
Nevertheless, as two of us pointed out in Ref. [3, the edge 
defects can be viewed as halfvortices and have nontriv- 
ial topological charge related to the winding number of 
vortices in the bulk. 

For a bounded flat nanomagnet, the winding number 
of vortices in the bulk is not a conserved quantity. This is 
illustrated by an example shown in Fig. [1] where a bulk 
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vortex with winding number n = ±1 is absorbed into the 
edge. Conservation of topological charges can be restored 
by assigning winding numbers to edge defects. In this 
case there are two such kinks at the edge of the film. The 
process shown in Fig.[T]then expresses the annihilation of 
a ±1 bulk vortex with two — i edge defects. Numerical 
simulations exhibiting similar annihilation of bulk vortex 
with edge defects can be found in Ref. [^. 

The winding number of a single edge defecj^ is defined 
as the line integral along the boundary dfl: |5| 

'-9r)-dr^±^. (1) 



27r Jen 



Examples of edge defects with half-integer winding num- 
bers are shown in Fig. [31 For a closed boundary the sum 
of the winding numbers of edge defects is also given by 
the above integral, but instead of integrating around one 
edge defect, the integral is carried out along the entire 
boundary. It was shown in Ref. |5| that this integral is 
related to the sum of winding numbers of vortices in the 
bulk. In general, for a film with g holes, we obtain 



edge 



bulk 

E' 



= 1 



(2) 



Here the winding numbers nt are integers for bulk de- 
fects and half-integers for edge defects. This conservation 
law has inrportant implications for the dynamics of mag- 
netization in nanomagnets [5|. Since defects with large 
winding numbers carry significant magnetic charge and 
are unfavored energetically in flat nanomagnets, most 
of the intricate textures observed involve only bulk vor- 
tices with winding number 71 = ±1 and edge defects with 

Topological considerations also place important con- 
straints on the possible structure of the domain walls in 
magnetic nanostrips JQl . Examples of such domain walls 
are shown in Figs. [3] and [S] Since edge defects are kinks 
of magnetization along the boundary, a domain wall in 
a magnetic strip must contain an odd number of kinks 
at each edge. Furthermore, the angle of magnetization 
rotation along the two edges must be compensated by 
the winding number of the bulk. Consequently, the total 
topological charge including contributions of vortices and 
edge defects must be zero in a head-to-head domain wall. 
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FIG. 1: A vortex (n = +1) absorbed by the edge can be 
viewed as its annihilation with two — i edge defects. The 
annihilation results in a uniform magnetization pointing to 
the right. 



The edge defects in nanomagnets are analogs of the 
so called boojums which exist at the surfaces and in- 
terfaces of superfluid "^He [SJ, l9J. In general, "boojum" 
refers to a topological defect that can live only on the 
surface of an ordered medium [10|. Boojums were also 
predicted and observed in some liquid crystals [ll|. An 
interesting system which is closely related to our study 
of flat nanomagnets is the smectic C films. Convention- 
ally, the ordered state of a liquid crystal is described by 
a unit vector n pointing along the long axis of the con- 
stituent molecules. Because order parameters n and — n 
are equivalent, it is known that vortices with half-integer 
winding numbers are allowed to exist in the bulk of ne- 
matic liquid crystals. On the other hand, for smectic C 
liquid crystals, the in-plane ordering of molecular orien- 
tations is described by an additional 2D unit vector c 
lying in the smectic planes and pointing to the tilt direc- 
tion of fi [l2l. Because rotating a tilted molecule by 180° 
around the normal of smectic layers does not return it 
to its original configuration, this unit vector c, like the 
magnetization ih, is a true vector. As discussed below, 
the vector nature of the order parameter rii is important 
to the confinement of halfvortices at the edge. 

In this article we review the structure and energetics of 
elementary topological defects in nanomagnets. In con- 
trast to the determination of the configuration of topo- 
logical defects in superfluids or liquid crystals where the 
energy is dominated by short range interactions, finding 
solutions of the vector field m(r) for topological defects 
in nanomagnets is considerably more difficult due to the 
nonlocal nature of dipolar interaction. We approached 
this problem from two opposite limits dominated by the 
exchange and dipolar interactions, respectively; the re- 
sults are presented in Sections lU and Hill Edge defects of 
smectic C films are discussed in Section lTVl where we also 
point out the similarities and differences of the two mod- 
els. We conclude with a summary of our major results in 
Section IVl 



II. EXCHANGE LIMIT OF FLAT 
NANOMAGNETS 

The magnetic energy of a ferromagnetic nanoparti- 
cle has two major contributions: the exchange energy 
A I |Vmp (ff'r and the dipolar energy (a<o/2) / |Hp d^r. 
The magnetic field H is related to the magnetization 
through Maxwell's equations, V x H = and V • (H + 
M) = 0. Here we disregard the energy of anisotropy, 
which is negligible for soft ferromagnets such as permal- 
loy. 

Analytical treatment of topological defects is gener- 
ally impossible due to the long range nature of dipolar 
interaction. One usually minimizes the energy numeri- 
cally to find stable structures of the magnetization field. 
Nevertheless exact solutions are possible in a thin-film 
limit [H, Hi]: i < w < A^/t < w log {wjt) defined for a 
strip of width w and thickness t. Here A = y/AjJioM^ is 
a characteristic length scale of exchange interaction. In 
this limit the magnetization only depends on the in-plane 
coordinates x and y, but not on z. More importantly, the 
magnetic energy becomes a local functional of magneti- 
zation [il,[i3: 



E[ih{r)]/At= / \Wm\^(fr+{l/A) 



f (m-fi)= 
Jan 



dr. (3) 



Here il is the two-dimensional region of the film, dfl is its 
line boundary, n _L r is unit vector pointing to the out- 
ward normal of the boundary, and A = 47rA^/tlog(w/t) 
is an effective magnetic length in the thin-film geometry. 
Eq. ^ is the familiar XY model [l| with anisotropy at 
the edge resulting from the dipolar interaction. Denot- 
ing m = (cos^, sin0), minimization of the energy E with 
respect to 9 yields the Laplace equation \/^9 = in the 
bulk and boundary condition n • V0 = sin 2 {9 + 9e) /A 
at the edge. 

Topological defects that are stable in the bulk are or- 
dinary vortices with integer winding numbers, which are 
well known in the XY model ^1]. The boundary term 
of model ^ introduces yet another class of topological 
defects that have a singular core outside the edge of the 
system. To be explicit, consider an infinite semiplane 
y > 0. Solutions satisfying the Laplace equation in the 
bulk and the boundary condition dy9 = sin26'/A at the 
edge y = are 0, ^^ 



tan6'(x,y) — ± 



y 



A 



x-X 



(4) 



The singular core is at {X, —A), distance A outside of the 
edge. Fig. [2]shows the magnetization fields of Eq. ((4]). As 
can be easily checked using Eq. ([T]) the winding numbers 
of these solutions are ±i, respectively. The halfvortices 
can not live in the bulk: as its singular core is moved 
inside the boundary, a string of misaligned spins occurs 
which extends from the core of halfvortex to the bound- 
ary [5|. The edge thus provides a linear confining poten- 
tial for the halfvortices. 



In the limit A/w — >• 0, magnetization at the edge is 
forced to be parallel to the edge, m — ±t. By exploiting 
the analogy between XY model and 2D electrostatics, 
one can use the method of images to deal with the ef- 
fects introduced by the boundary [ll Q • In this analogy 
the vortex is mapped to a point charge whose strength 
is given by the corresponding winding number. However, 
unlike the electrostatics, the "image" charge induced by 
the boundary has the same sign as the original. The 
above solution ^ with A = looks just like a n = ±1 
vortex with its core sitting at the edge. The assignment 
of half- integer winding number n = ±i to the edge de- 
fect is thus consistent with the electrostatics analogy in 
the sense that the winding number is doubled by the re- 
flection at the edge [5'|. 

An exact solution for a domain wall was also obtained 
in this limit 5]. Consider a strip \y\ < w/2. It has 
two ground states with uniform magnetization: 9 — or 
TT. Domain walls interpolating between the two ground 
states are given by 



tan9{x,y) — ± 



cos ky 



sinhfc(x — X) ' 



(5) 



where the wavenumber k ^ t:/{w + 2A). The magneti- 
zation field of Eq. ([5]) (shown in Fig. |3]) is reminiscent of 
the so called 'transverse' domain walls (Bottom panel of 
Fig. [3]) observed in micromagnetic simulations [l5|. 

Unlike domain walls (kinks) in Ising magnet, the do- 
main wall described by Eq. ([5]) is a composite object con- 
taining two edge defects with opposite winding numbers 
±^. The singular cores of the two halfvortices reside 
outside the film, a distance A away from the edges. One 
can understand the stability of the domain wall using 
the electrostatics analogy: the attractive 'Coulomb' force 
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2: Edge defects with winding numbers n = -|-i (top) 
-| (bottom) in the exchange limit. 
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FIG. 3: Top: Magnetization of the head-to-head domain wall 
solution (O . It is composed of two edge defects with opposite 
winding numbers ±i. Bottom: a transverse domain wall ob- 
served in a micromagnetic simulation using OOMMF [20| in 
a permalloy strip of width w = 80 nm and thickness f = 20 



pulling together the two halfvortices is balanced by the 
confining force from the edges. 

The total energy of the domain wall solution Eq. ([S]) 
evaluates to E k, 27r^i(l -I- log(u'/7rA)). As expected 
for the XY model, the exchange energy depends loga- 
rithmically on the system size which is the width of the 
strip vj in our case. It also depends logarithmically on a 
short distance cutoff which is provided by A here. After 
restoring the energy units and expressing A in terms of 
the relevant parameters, we obtain the following domain 
wall energy in the exchange limit 



^ „ . , /ewilog(w/t) 

E^DW-STT^ilogf ^^ ' ' 



Try 



(6) 



The energy depends linearly on the thickness of the film 
t and only weakly (logarithmically) on the width. These 
relations are important to the understanding of the hys- 
teresis curves of asymmetric magnetic nanorings [l6| . 



III. DIPOLAR LIMIT OF FLAT 

NANOMAGNETS 

The thin-film limit discussed in the previous section 
is inaccessible to most experimental realizations of nano- 
magnets, in which the dipolar interaction is the primary 
driving force. In this section we discuss the structure 
and energetics of topological defects and domain walls 
in the opposite limit where the energy is dominated by 
the dipolar interaction. Our strategy here is first to 
find structures which minimize the magnetostatic energy 
{fio/2) J jHpd'^r and then to include exchange interac- 
tion as a perturbation. However, energy minimization 
in the dipolar limit is relatively difficult due to following 
reasons. Firstly, as opposed to the local exchange inter- 
action, the dipolar interaction is long-ranged. Secondly, 
in many cases the magnetostatic energy has a large num- 
ber of absolute minima. One thus has to search among 
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FIG. 4: an antivortex (top), and a — ^ edge defect (bottom) 
in the dipolar limit. 



these minima for one with the lowest exchange energy, 
making it a degenerate perturbation problem. 

The magnetostatic energy of a given magnetization 
field m(r) can also be expressed as the Coulomb interac- 
tion of magnetic charges with density p^ (r) = — AfoV-rh, 
where Mq is the saturation magnetization. Being positive 
definite, the magnetostatic energy has an absolute min- 
imum of zero, which corresponds to a complete absence 
of magnetic charges. A general method to obtain the ab- 
solute minima of magnetostatic energy was provided by 
van den Berg in 1986 [l7|. For magnetic films with arbi- 
trary shapes, his method yields domains of slowly varying 
magnetization separated by discontinuous Neel walls. In 
the following we look for structures that have the desired 
winding number and are free of magnetic charges, i.e. 
V • m = in the bulk and n • rii = on the boundary. 

We start by examining the vortex solutions of XY 
model. In polar coordinate, a vortex with winding num- 
ber n is described by 6{x,y) — ncj) + 9q, where 9o is 
a constant and (j) — arctan(y/a;) is the azimuthal an- 
gle. Among these solutions, only the n = 1 vortex with 
^0 — ""/^ has zero charge density and survives in the 
dipolar limit. Its energy then comes entirely from the ex- 
change interaction and diverges logrithmically with sys- 
tem size R: E K, 2'KAt\og{R/\). Here the short distance 
cutoff is given by the exchange length A. 

The antivortex solutions of the XY model always carry 
a finite density of magnetic charge and thus are not a 
good starting point to obtain the n = —1 defect in the 
dipolar limit. Fortunately, a magnetization field with 
winding number —1 and free of bulk charges is realized 
by a configuration known as the cross tie [y, |l^ (top 



panel of Fig. |4]). It consists of two 90° Neel walls normal 
to each other and intersecting at the singular core. The 
magnetization field of an antihalfvortex (winding num- 
ber — ^) is obtained by placing the core of a cross tie at 
the edge of the film (bottom panel of Fig. |4]). Since the 
magnetization along the edge is parallel to the boundary, 
the structure is also free of surface charge. As one moves 
from left to right along the edge the magnetization ro- 
tates counterclockwise through tt. This is in agreement 
with the definition ([1]) for an antihalfvortex. 

The energy of an antivortex or an antihalfvortex grows 
linearly with the length of the Neel walls L emanating 
from it: 



E ^ (jtL + Er. 



(7) 



The surface tension of the wall a has contributions from 
both exchange and dipolar interactions. In magnetic 
films with thickness exceeding the Neel-wall width (of 
order A), it is given by [6| 



(8) 



(7 = 2V2(sin 6*0 - 6*0 cos 6*0) A/X, 



where 29q is the angle of magnetization rotation across 
the wall. In thinner films (i ^ A) the magnetostatic 
term becomes substantially nonlocal and the Neel walls 
acquire long tails [l9| . 

There is no charge- free configuration for the +^ edge 
defect. In addition, one could also observe from micro- 
magnetic simulations that most of the magnetic charges 
of a transverse domain wall in a strip is accumulated 
around the +^ defect. Thus, in the dipolar limit, the 
chargeful +i defect is prone to decay into a — ^ edge 
defect and -|-1 vortex in the bulk. We next turn to the 
discussion of the structure of domain walls in this limit. 

An intrinsic problem arises when one tries to apply van 
den Berg's method to find the structure of domain walls. 
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FIG. 5: Top: a magnetization configuration free of bulk mag- 
netic charges, —V • M = 0, and containing two —1/2 edge 
defects and a -f 1 vortex in the middle. Parabolic segments 
of Neel walls are shown by dashed lines. Bottom: a head- 
to-head vortex wall obtained in a micromagnetic simulation 
using OOMMF [2C| in a permalloy strip of width w — 500 nm 
and thickness i = 20 nm. 



That is because a head-to-head domain wall carries a 
fixed nonzero amount of magnetic charge: Qm = 2MQtw. 
However, these magnetic charges tend to repel each other 
and spread over the surface of the sample, much the 
same as the electric charges do in a metal. Based on 
this principle, we provided in Ref. 7] a construction of 
the head-to-head domain wall that is free of bulk mag- 
netic charges. All of the charge Qm is expelled to the 
edges. The resulting structure is shown in the top panel 
of Fig. [SI It resembles the structure known as the 'vor- 
tex' domain wall (bottom panel of Fig. [5|) predicted to be 
stable in regimes dominated by dipolar interaction |15| . 
Both structures contain two — ^ edge defects sharing one 
of their Neel walls and a -1-1 vortex residing at the mid- 
point of the common wall. 

The variational construction contains charge-free do- 
mains with uniform and curling magnetization separated 
by straight and parabolic Neel walls. In a strip |y| < w/2, 
the two — i edge defects share a Neel wall x — y where 
the vortex core {v,v) is located. The two curling do- 
mains in the regions ±i; < ±y < w/2 are separated by 
parabolic Neel walls {x — v)'^ = {2y ± 'w){2v ± w) from 
domains with horizontal magnetization; they also merge 
seamlessly with other uniform domains along the lines 
X = V and y — v. 

The location (w, v) of the -1-1 vortex on the shared Neel 
wall is a free parameter of our variational construction. 
The structure remains free of bulk charge as the vortex 
core moves along the diagonal x = y. When it reaches 
one of the edge, its annihilation with the — ^ edge defect 
creates a widely extended +-^ edge defect. The resulting 
structure (top panel of Fig. |n| is topologically equivalent 
to the transverse domain wall (bottom panel of Figs. [3] 
and [6]) discussed in the previous section. 

The equilibrium structure for given strip width w and 
thickness t is determined by minimizing the total energy 
of the composite domain wall with respect to vortex coor- 
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FIG. 6: Top: A model vortex wall when the vortex is absorbed 
by the edge forming an extended -\-\ edge defect along the 
upper boundary. Bottom: Transverse wall observed in micro- 
magnetic simulation. 



/'""" 
/ 


\ 






/' 
/ 
/ 

/ 
y' ^-'' 


\ 




" -^ ^ ~ — - ll- 


^1^ 


-^'— ' 


'L. ^- -^ ^ ^ ^ 




/ 


x^ 






^^^ 


\ 

\ 


/ 








t=20 nm - 


_\ 


t=5 nm — 


-- \ 







1.3 
1.2 \ 
1.1 
S 1 

LU 

2" 0.9 

0.8 

0.7 

0.6 

-0.5 -0.4 -0.3 -0.2 -0.1 0.1 0.2 0.3 0.4 0.5 

V 

FIG. 7: Energy of the vortex domain wall as a function of the 
vortex position v aX & fixed strip width ui = 50 nm for several 
thicknesses t. 



dinate v. The total energy contains the following terms, 
(a) The exchange energy of the two curling domains fi 
around the vortex core. It is given by At J^{W9)^(Pr 
and of the order At\og{w/X). (b) The energy of Neel 
walls, which can be computed as a line integral t J cr{i) d£, 
where d£ is a line element of the wall. The surface tension 
a, which depends on the angle of spin rotation across the 
wall, is given by Eq. (|8]). This term is of order Atw/X. 
(c) The magnetostatic energy coming from the Coulomb 
interaction of magnetic charges spreading along the two 
edges. It is of the order Aw{t'^/X^) log(w/i). 

By combining the above three contributions, the to- 
tal energy curve E(v) for a fixed width w and varying 
thickness t is shown in Fig. [T) For substantially wide 
and thick strips, the curve attains its absolute minimum 
as the vortex is in the middle of the strip, in agreement 
with numerical simulations [15|. A local minimum de- 
velops with the vortex core at the edge of the strip as 
the thickness decreases. This solution corresponds to the 
transverse wall shown in Fig. [S] The transverse wall be- 
comes the absolute minimum as the thickness is further 
reduced and the vortex wall {v = 0) is locally unstable. 
It should be noted that the above calculation for thin 
films, e.g. t = 1 nm, is only an extrapolation. For films 
with small cross section (but not in the exchange limit), 
our variational approach can not be trusted. Nonethe- 
less, the method is illustrative and indeed shows that the 
three-defects wall structure is unstable when approaching 
the exchange limit. 



IV. HALFVORTICES IN SMECTIC FILMS 

The XY model, applicable in the thin- film limit, pre- 
serves a symmetry between topological defects with op- 
posite winding numbers, namely, ±1 vortices have ex- 
actly the same energy in model (|31) (so do ±i edge de- 
fects) . Since vortices of opposite winding numbers carry 



different magnetic charges, the degeneracy is lifted in 
thicker and wider strips where the dipolar interaction 
becomes more important. The configuration of the topo- 
logical defects in the extreme dipolar limit discussed pre- 
viously clearly shows this asymmetry. One can also break 
this symmetry by assigning different penalties to splay 
(V • rii 7^ 0) and bending (V x m 7^ 0) deformations: 

E[m{x)] = I [Ki{\7 -xaf + K2{^ y.mf]Sr 
Jn 

+ (1/A) / (n-m)2dr, (9) 

Jdn 

Here the elastic constants Ki and K2 have energy unit, 
whereas the edge anisotropy 1/A scales as the inverse 
length times energy. With the unit vector m identified as 
the c-director field, this energy functional also describes 
the elastic energ y of a chiral smectic film f2r| or a Lang- 
muir monolayer (22| with planar boundary conditions. 

The case Ki = K2 corresponds to the XY model and 
the exchange limit discussed in Sec.|TTl By choosing Ki > 
K2 we discourage splay, which is similar to a penalty 
for magnetic charges in the bulk. The dipolar limit thus 
corresponds to the regime where the bend energy is small 
compared to those of splay and edge anisotropy. In what 
follows we focus on the extreme dipolar limit K2 = 0. 

First, the -1-1 vortex solution 9{r) = (f) + tt/2 remains 
an energy minimum of model © for arbitrary Ki and 
K2. The -hl/2 edge defect in the XY limit, Eq. @ with 
the '-I-' sign, also is a stable configuration for arbitrary 
K and A except that the singular core is pushed further 
outside the boundary, a distance (1 -I- e)A away from the 
edge. Here e = (i^i - i^2)/(i^i + K2). 

Since the bulk term of the energy functional (JH) does 
not have an intrinsic length scale, an exact scale-invariant 
solution for an antivortex has been obtained in the dipo- 
lar limit K2 = 0: 



9{x,y) = (j)- 



1(72: 



(10) 



where cj) — arctan(y/a;) is the azimuthal angle. This so- 
lution is singular aX 4> = ±7r/4, where the first derivative 
d9/d(j) diverges. A complete solution of the antivortex 
nevertheless can be obtained by continuing the above so- 
lution outside of the interval \(j)\ < tt/A periodically. The 
result is shown in Fig. |Sl^a). 

Analytical solutions of antihalfvortex for arbitrary A 
are yet to be found. In the limit A — ^ achieved in 
boundaries with very strong anchoring force, the unit 
vector m is forced to be parallel to the edge. In this 
limit, the — i edge defect can be constructed following 
the same trick for the antihalfvortex in the dipolar limit 
of nanomagnets. The resulting configuration is shown in 
Fig. IHl^b). Compared with their counterparts in the XY 
model, the antivortex and antihalfvortex in Fig. [5] are 
closer to the cross tie configuration (or half of it) shown 
in Fig. m 

Although topological defects of the generalized elastic 
model show some similarities with those of the magnetic 
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FIG. 8: an antivortex (top), and a — ^ edge defect (bottom) 
in the elastic model with K\ = 1, K2 ~ 0, and A = 0. 



problem in the dipolar limit, there is an important dif- 
ference regarding the scaling of their energy with sys- 
tem size. Since solution (fTO|) is scale-invariant, the en- 
ergy of the defects also diverges logarithmically: E ~ 
const X A'l log (i?/a). Here a is a short distance cutoff. 
In the case of — ^ edge defect, a is of the order of A. 
However, as discussed in the previous section, the energy 
of both antivortex and antihalfvortex scales linearly with 
the length of Neel wall. The nonlocal dipolar interaction 
in the magnetic problem results in a natural length scale 
A = ^A/fioAP. By contrast, there is no such length 
scale in the elastic model ©, so the dependence of the 
energy on the system size is logarithmic. 



V. CONCLUSION 

We have discussed the topological properties of edge 
defects in XY systems with a broken 0(2) symmetry 
at the boundary. In particular, we discussed two physi- 
cal systems containing such edge defects: the 2D smec- 
tic C films and nanomagnets with a planar geometry. 
Since spins at the boundary have two degenerate pre- 
ferred directions, i.e. parallel or antiparallel to the tan- 
gent of boundary, the edge defects are manifested as kinks 
of magnetization along the edge. Moreover, they carry 
half-integer winding numbers and thus can be viewed as 
halfvortices confined to the edge. Conservation of topo- 
logical winding number can only be established by in- 
cluding contributions from the edge defects. As we have 



pointed out before [3431 ; edge defects should be included 
along with ordinary vortices as the elementary topolog- 
ical defects in flat nanomagnets. Indeed, domain walls 
which play an significant role in the dynamics of magnetic 
nanostrips and nanorings are composite objects consist- 
ing of two or more of these elementary defects. 

Analytical solutions of halfvortices and transverse do- 
main walls were obtained in a thin-film limit where the 
exchange interaction is the dominant force determining 
the shape of topological defects. The magnetic problem 
is reduced to the familiar XY model with an anisotropy 
at the edge. Domain walls stable in this regime are com- 
posed of two edge defects with winding numbers ±i. By 
analogy with 2D electrostatics, the stability of transverse 
domain wall can be understood as resulting from a bal- 
ance of the attractive Coulomb force between the oppo- 
sitely charged halfvortices and the confining force from 
the edges. 

Energy minimization is relatively difficult in the oppo- 
site limit dominated by the nonlocal dipolar interaction. 
Nevertheless, by focusing on structures which are free of 
bulk magnetic charge, we are able to find structures of 
topological defects stable in this regime. The -1-1 vortex 
of XY model with circulating magnetization remains a 
stable defect in the dipolar limit. The —1 vortex sur- 
vives in this limit but is severely deformed; it has the 
cross tie structure consisting of two 90° Neel walls in- 



tersecting at the singular core. The configuration of the 
— i edge defect is constructed by placing the core of a 
cross tie at the boundary. The +^ defect carries a finite 
amount of magnetic charge and is unstable in this limit. 

We have presented a variational construction of the 
vortex domain wall which is composed of two — ^ edge 
defects and a -t-1 vortex. By varying the location of the 
center -t-1 vortex, the construction interpolates between 
the vortex wall and the transverse wall. Variational cal- 
culation of the domain wall energy reveals that the vor- 
tex wall is indeed stable in the dipolar limit whereas it 
becomes an energy maximum in thin and narrow strips. 

Finally, we have discussed structures of topological 
defects in an elastic model which generalizes the XY 
model of the thin-film limit. Calculations in this model 
are simplified by the replacement of non-local interac- 
tions between magnetic charges by a term that penalizes 
the existence of magnetic charge in a local fashion. This 
model is applicable to smectic C films, but may provide 
insight into magnetic configurations. In particular, the 
allowed topological defects are the same in both systems. 
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